Draft version February 5, 2008 

Preprint typeset using I^T^X style cmulateapj v. 10/09/06 



O 
O 

(N 

u 

< 
in 

(N 

> 
oo 

OS 

in 
m 
O 

o 

6 ■ 
a: 



13 



SPECTRAL EVOLUTION OF NGC 1313 X-2: EVIDENCE AGAINST THE COOL DISK MODEL 
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ABSTRACT 

The presence of a cool multicolor disk component with an inner disk temperature kT = 0.1 ~ 
0.3 keV at a luminosity L > 10 40 ergs s _1 has been interpreted as evidence that the ultraluminous 
X-ray source NGC 1313 X-2 harbors an intermediate-mass black hole (IMBH). The temperature of 
a disk component should vary with luminosity as L oc T 4 . However, upon investigating the spectral 
evolution with multiple XMM-Newton observations, we found that the cool disk component failed 
to follow this relation with a confidence level of 0.999964. Indeed, the luminosity decreases as the 
temperature increases, and the luminosities at high temperatures are more than an order of magnitude 
less than expected from the L oc T 4 extrapolation of luminosities at low temperatures. This places 
a strong constraint against the validity of modeling the X-ray spectra of NGC 1313 X-2 as emission 
from the accretion disk of an IMBH. The decrease in luminosity with increasing temperature of the 
soft component follows the trend suggested by a model in which the soft emission arises from an 
outflow from a stellar-mass black hole with super-Eddington accretion viewed along the symmetry 
axis. Alternatively, the spectra can be adequately fitted by a p-free disk model with kT m 2 keV 
and p f» 0.5. The spectral evolution is consistent with the L oc T 4 relation and appears to be a high 
luminosity extension of the L — kT relation of Galactic black holes. This, again, would suggest that 
the emission is from a super-Eddington accreting stellar mass black hole. 

Subject headings: black hole physics - accretion, accretion disks - X-rays: binaries - X-rays: galaxies 
- X-rays: individual (NGC 1313 X-2) 



1. INTRODUCTION 

Ultraluminous X-ray sources (ULXs) are pointlike, 
nonnuclear X-ray sources with luminosities above the 
Eddington limit of a 20 M m black hole (3 x 10 39 ergs 
s ; for a review, see iFabbiano fe Whitdl200l . Strong 
variability from many ULXs indicates they are compact 
objects. If the X-rays are emitted isotropically below 
the Eddington limit, then the hig h luminosity would in- 
dicate that ULXs contain IMBHs ([Colbert fe Mushotzkvl 
Il999t [Kaaret e t al. 2001). However, the emi ssion could 
be mechanically or relat ivistically beamed (|King et al.l 
l200lHKording et alj2002ft. or exceed th e Eddington limit 
([Begelmanl l2002t lEbisawa et al.l l2003[ k in these cases, 
IMBHs are not required. 

The sum of a power-law and a multicolor disk (MCD; 
iShakura fc Sunvaevl 119731 ) has been found to describe 
well the X-ray spectra of black hole binaries. Fitting 
ULX spectra with the sam e model often show s a disk 
component, first reported in lKaaret et alj (|2003t ). with a 
disk temperature around 0.1-0.4 keV which is much lower 
than that found in stellar-mass black holes. The cool- 
disk high-luminosity regime sugges ts that ULXs har- 
bor IMBHs with masses of 10 3 M w (iKaaret et all 120031 : 
iMiller et al.ll2003t iFenTfc Kaared[20rM ~ 

The radial temperature of an MCD follows a power- 
law form as kT oc R~ p , where p = 0.75. When tak- 
ing into account advective cooling, the po wer-law depen- 
dence p is predicted to decrease until 0.5 ([Watarai et al.l 
2000), which usually happens at very high accretion 
rates when the source enters the 'slim disk' or opti- 
cally thick advection-do minated accretion disk phase 
(|Abramowicz et al"1ll988| ). The p-free model is similar to 
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the M CD model but allows p to vary (Mincshig e et al.l 
1994), and has been successfully applied in inter- 
preting the very high s tate of black hole binarie s 
such as XTE J1550-56 4 ([Kubota fc Makishimal[200l . 
IVierdavanti et all (2006) found that four ULXs that were 
previously reported to contain cool disks, could also be 
adequately fitted using the p-free model with disk tem- 
peratures of around 2 keV and p « 0.5. They appear in 
the L — kT diagram as an extension of stellar-mass black 
holes with super-Eddington radiation. 

NGC 1313 X-2 is one of the UL Xs with a spectrum 
that can be modeled by a cool MC D ([Miller et alJl2003f >. 
a hot MCD (IStobbart et alj|2006ft. or a hot p-free disk 
([Mizuno et al.l l2007t IVierdavanti et~aT1 l2006f ) . The to- 
tally opposite physical consequences derived from the 
cool and hot disk models inspire us to investigate whether 
one can be ruled out. The source resides in an opti- 
cal nebular supershell (Pakull fc MirionH |2003[ ) , whose 
kinetic energy is much larger than a typical supernova 
explosion. The shell is either a hypernova remnant or a 
continually powered nebula possibly driven by the out- 
flow from the ULX. 

In this Letter, we revisited t he 12 archival XMM- 
Newton observations described in lFeng fc Kaaret] ([2006) . 
We present the spectral analysis in Sj2]and discuss the na- 
ture of the ULX in fl3] We adopt a dist ance to the host 
galaxy of 4.13 Mpc ([Mendez et ai1l2002D . 

2. SPECTRAL ANALYSIS 

iFeng fc Kaaretl (|2006h analyzed 12 XMM-Newton ob- 
servations of NGC 1313 and described the spectral state 
transitions of the two ULXs X-l and X-2. This previ- 
ous Letter identified a tight, positive correlation between 
the luminosity and the power-law photon index of X-l, 
while X-2 presented an opposite pattern of high/hard 
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Fig. 1. — 1 cr confidence contours between the absorption col- 
umn density TVh and the inner disk temperature kT for the MCD 
m odel. Difference colors represent for different observations noted 
in IFeng & Kaaret (2006). The plus signs indicate the best-fit val- 
ues. For each observation, the range of kT in a contour is adopted 
as its error, used in Fig. [2] 
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Fig. 2.— L - kT diagram of NGC 1313 X-2 for the MCD model, 
with total 0.3—10 keV luminosity (top) and 0.3—10 keV disk lumi- 
nosity (bottom). Solid lines indicate the best-fit L oc T 4 relation, 
with x 2 /dof = 28/5 for the L — kT dia gram and x 2 /dof = 41/5 for 
the Ldisk — kT diagram. Dashed lines are the best-fit L oc T n rela- 
tion, where n = —3.7 ±0.7 for the L — kT diagram and n = —8 ±2 
for the Ldisk — kT diagram. 

to low/soft spectral evolution consistent with previous 
ASCA observations (jZampieri et al.ll2004T) . 

In 6 of the 12 observations, NGC 1313 X-2 exhibits a 
significant soft excess which can be fitted by an MCD 
model with temperatures of 0.1-0.3 keV, contributing 
20%-60% to the total 0.3-10 ke V flux. All spec t ral pa - 
rameters are listed in Table 1 of IFeng fc Kaaretl ([2006) . 
Because the MCD component peaks where the interstel- 
lar absorption is of high importance and the observed 
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Fig. 3.— L - kT diagram of NGC 1313 X-2 for the blackbody 
model. The dashed line indicates the best-fit L oc T n relation, 
where n = -3.1 ± 0.5 and x 2 /dof = 2.7/4. 

MCD flux is often low relative to the power-law flux, we 
have carefully reevaluated the error on the temperature 
kT to include the influence of the absorption column den- 
sity ./Vh, rather than usi ng the default 1 a valu e derived 
from the spectral fits in IFeng fc Kaaretl (|2006l ) . We cre- 
ated two-dimensional plots of the 1 a confidence contours 
allowing both Ng and kT to vary (see Fig. [T]). The con- 
tours are calculated on a grid of points with step sizes 
of 0.005 keV and 5 x 10 19 cm~ 2 . To find the error on 
kT, we use the full range in kT covered by the 1 a confi- 
dence range, i.e. from the leftmost to the rightmost point 
of the contour. This error includes the influence of JVh 
and is larger than the conventional definition of the error 
derived from the error command in XSPEC. For these 
six observations, the 0.3-10 keV total luminosity versus 
temperature is plotted in Figure[2](top), while the 0.3-10 
keV disk luminosity versus temperature is plotted in the 
bottom panel. 

We refitted all of the observations using a comptt 
model in place of the power-law, and the MCD tempera- 
ture changed by less than 0.02 keV in all cases. Replac- 
ing the power-law model by a Comptonization model has 
little effect on the MCD component, because the source 
spectrum does not exhibit cur vature or a break at high 
energies (jStobbart et al.l 120061 ). We also fitted the data 
with a model consisting of a simple blackbody and a 
power-law. The results are shown in Figure [3] The er- 
ror on the temperature is derived in the same way as 
for the MCD component, including the influence of the 
absorption on the error estimate. 

To see if the data are consistent with a L oc T 4 relation, 
we present a best- fit L cx T 4 line on each plot. Since the 
uncertainty on the temperature is much larger than on 
the luminosity, we treated L as the independent variable 
and T as the dependent variable and took into account 
the error on the temperature in the fitting. The best- 
fit L cx T A line results in a % 2 /dof = 28/5 (probability 
= 3.6 x 10~ 5 ) for the L—kT diagram and a x 2 /dof = 41 /5 
(probability = 9.4 x 10~ 8 ) for the Ldisk — kT diagram, 
respectively. The L cx T 4 relation is strongly excluded 
by the data in both cases. We note that if we use the 
la error derived from the error command in XSPEC, in 
which the Ng. influence is not considered, the the best-fit 
L cx T 4 relation will result in a x 2 /dof = 68/5 (proba- 
bility = 3 x 10" 13 ) for the L - kT diagram. 

We also performed a fit where the power-law index 
among L and T was left as a free parameter, n, with L oc 
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T n . The best- fit L oc T n relations are shown as dashed 
lines in Fig.[2]with n = — 3.5±1.0 for the L — kT diagram 
and n = — 7±3 for the L^isk — kT diagram. These provide 
reasonable fits, with x 2 /dof = 2.7/4 (probability = 0.61) 
for the L — kT diagram and x 2 /dof = 5.9/4 (probability 
= 0.21) for the Ldisk — diagram. 

The best-fit disk inner radius R ln varies a lot between 
observations with different temperatures. Assuming a 
face-on disk at a distance of 4.13 Mpc, the best-fit i?; n 
is about (2 — 4) x 10 4 km for the low temperature group 
(kT = 0.11 - 0.12 keV), but around (1 - 2) x 10 3 km for 
the high temperature group (kT = 0.22 - 0.24 keV). 

We fitted all spectra of X-2 with the p-free model. Due 
to the limited number of photons at high energies, errors 
on the temperature from a single observation are large, 
in particular for observations in the low luminosity state. 
We combined the 12 observations into 4 groups on basis 
of similar spectral properties, and performed simultane- 
ous spectral fitting for each group to reduce the errors. 
Best- fit parameters for the 4 groups are listed in Table. [TJ 
We plotted the L — kT diagram for p-free model in Fig- 
ure |U The best-fit L oc T 4 line is shown as a solid line 
in Figure [U Although the 4 data points are insufficient 
to obtain an unambiguous pattern of the spectral evolu- 
tion, they are roughly consistent with aiocT 4 relation, 
X 2 /dof = 9/3. In Figured] (bottom), we overplotted the 
L — kT patterns of three Galactic black hole binaries: 
GX 339-4, LMC X-3, an d XTE J1550-564 (data from 
iGierlihski fc T5o^l2004ah . 

We also applied the MCD plus power-law model for the 
four grouped observations and found similar results with 
ungrouped data. Two of the four observations, group 1 
and 2, present significant MCD components with tem- 
peratures of 0.22 and 0.11 keV, respectively, correspond- 
ing to a total 0.3-10 keV luminosity of 3.6 x 10 39 and 
2.1 x 10 40 ergs s" 1 . A gain, the correlation between lumi- 
nosity and temperature is opposite that expected if the 
soft thermal component arises from an accretion disk. 

3. DISCUSSION 

The spectral evolution of NGC 1313 X-2 is inconsistent 
with the L oc T 4 relation expected if the cool thermal 
component is, indeed, emission from an accretion disk 
at a high level of significance. Detection of the cool disk 
component in ULXs is one of the arguments that support 
the presence of IMBHs. However, whether we can model 
the soft excess in the X-ray spectra of ULXs with a cool 
MCD, and even t he reality of the so ft exce ss itself, has 
been questioned. iGoncalves fc Sorial (l2006h found that 
in some bright ULXs the presence of the soft excess de- 
pended on how the power-law continuum was modeled. 
By selecting different energy ranges to fit the power-law 
continuum, the soft excess could be replaced by a soft 
deficit. Similar soft excesses have been found in man y 
active galactic nuclei ( AGNs; IGierlihski fc Do ne 2004b) , 
which can be phenomenologically modeled as blackbody 
emission with a temperature of 0.1-0.3 keV over a wide 
range of AGN luminosity and black hole mass. The disk 
emission in AGNs peaks in the ultraviolet band, unable 
to produce the soft excess peaked in the X-ray band. 
The soft excesses in ULXs are similar to those in AGNs, 
which suggests that the y could have the same origin 
(|Goncalves fc Sorial[2006l ). iKind (|2004D pointed out that 
mass outflow should occur in Eddington-limited accret- 
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Fig. 4.— Top: L - kT diagram of NGC 1313 X-2 for the 
p-free model from the four grouped XMM-Newton observations 
(for parameters, see Table [TJ. The solid line indicates the best-fit 
L oc T 4 relation. Bottom: Same as the top panel but with L — kT 
patterns of three Galactic black hole binari es: GX 339-4, LM C X-3, 
and XTE J1550-564 (data from IGierlihski fc Donj|2004al) . 

ing compact objects, and could produce the blackbody- 
like soft excess in ULXs and AGNs via a Compton thick 
envelope around the central object. Therefore, it is im- 
prudent to derive the mass of the compact object before 
we can determine the reality of the soft excess. 

Spectral evolution offers a new means to test whether 
the soft excess is truly disk emission. The MCD dom- 
inates the high/soft spectral state of black hole bina- 
ries. Its evolution follows aiocT 4 pattern robustly and 
has been observed in man y Galactic black hole binaries 
(jGierlihski fc Pond [2004ah . In Figure EJ either with the 
total luminosity or with the disk luminosity, the obvi- 
ous departure from the L oc T A relation with the MCD 
model places a strong constraint that the MCD model 
is not the correct interpretation of the soft excess in the 
spectra of NGC 1313 X-2. The expected L oc T 4 relation 
is very strongly inconsistent with the data. Luminosi- 
ties at high temperatures are more than one order of 
magnitude less than expected from the best-fit L oc T 4 . 
Although small deviations between L and T 4 have been 
observed in Galactic black h ole binaries, such large de - 
viations have not been seen (jGierlihski fc Dondl2004ah . 
The error on the temperature must be treated carefully 
because of the low flux of the MCD component and the 
fact that the MCD peak occurs where there is strong 
curvature in the spectrum due to interstellar absorption. 
Our quoted uncertainties in the temperature take into ac- 
count the possible effects of the absorption column den- 
sity. Even with these conservatively evaluated errors, the 
predicted L — T 4 relation is ruled out by the data at a 
confidence level of 4.1er. 
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TABLE 1 

Best-fit spectral parameters of NGC 1313 X-2 with the p-free model 



group 
(1) 


observations 
(2) 


N H 
(10 21 cm" 2 ) 
(3) 


kT 
(keV) 
(4) 


V 

(5) 


0.3-10 keV L 
(10 39 ergs s" 1 ) 
(6) 


x 2 /dof 
(7) 


1 


1, 8, 10, 11 


1.95±0.03 


2.0±0.1 


5000+ 0003 
u.ouuu_ 0000 


2.7±0.2 


986/838 


2 


3, 4, 9, 12 


2.18±0.07 


2.5±0.1 


n c;8Q+0-009 
O-t>Sy_ .004 


9.5±0.5 


2432/2420 


3 


6, 7 


2.22±0.05 


1.7±0.1 


500+ 002 
u.ouu_ 000 


3.6±0.3 


477/461 


1 


2, 5 


2.39±0.09 


2.5±0.4 


n 5m +0.006 

u.oui_ 001 


5.2±1.0 


560/409 



Note. — Col. (1): Group index for combined spectral fitting. Col. (2): Observation index quoted from Feng & Kaaret (2006). Col. 
(3): Column density along the line of sight. Col. (4): Disk inner temperature. Col. (5): Power-law dependence of temperature on radius 
with a lower bounce of 0.5. Col. (6): 0.3—10 keV intrinsic luminosity. Col. (7): \ 2 an d degree of freedom. All errors are at the 1 a level. 



The contraction of the disk inner radius at high 
temperatures is reminiscent of the "anom alous regime" 
in G alactic black holes GRO J1 655— 40 dKubota et al 
200l|) and XTE J1550-564 (|Kubota fc Makishim 



2004f ). where strong disk Comptonization takes place. 
However, the L oc T 4 correlation for these sources is 
maintained if the total luminosity is considered instead 
of the disk luminosity. NGC 1313 X-2 fails to follow the 
L oc T 4 for either the total o r the disk luminosity. 

iFabrika fc Mcsch ervakovl (|2001h have proposed that 
ULXs are similar to the Galactic Galactic X-ray binary 
SS 433, which has a powerful outflow with a kinetic 
power in the range 3 x 10 38 to 1 x 10 40 ergs s _1 , but 



viewed along the jet axis. iBegelman et alj (|2006f ) sug- 
gest that SS 433 is highly super-Eddington, by a factor 
of 5000, and that ULX s are also highly super-Eddington. 
iPoutanen et alj (|2007f) considered the X-ray spectra in a 
similar model, also motivated by analogy with SS 433, 
and proposed that ULXs are geometrically beamed by 
a factor of 2-7, the soft excess is due to thermal emis- 
sion from the sphcrization radius where the outflow orig- 
inates, and the harder emission is from the accretion 
disk. Quasi-simultaneous observations made by XMM- 
Newton (part of these observations) and optical tele- 
scopes in 2004 showed that the optical counterpart of 
NGC 1313 X-2 did not brigh ten during the X-ray out- 
burst (jMucciarelli et aLll2006h . which might suggest that 
the outer region of the disk does not see the X-ray emis- 
sion and that the X-rays are geometrically beamed. 

In this scenario, the luminosity is not proportional to 
the temperature of the soft excess with a power of 4 as 
in the sub-critical accretion regime, instead, it varies in- 
versely with the temperature of the soft excess, i.e., the 
temperature at the spherization radius T sp . Our spec- 
tral results for N GC 1313 X-2 are cons istent with the 
outflow model bv IPoutanen et a l. (2007), which show a 
correlation between L and T with a negative power-law 
index (Figure^. Under this scenario, the soft excess is 
supposed to be modeled by a blackbody model instead 
of an MCD model. We thus fitted spectra for the six ob- 
servations with an absorbed blackbody plus power-law 
model, which also showed adequate fits. The L — kT dia- 
gram for the blackbody model is shown in Figure [3] The 
blackbody temperatures of around 0.1-0.2 keV indicate 
an accretion rate of about 20-100 times the Eddington 



rate and a spherization radius of about 200-1000 GM/c 2 . 
The best-fit dashed line indicates a power-law index of 
—3.2 between the two groups of data. The trend between 
L and T is as expected for a super-Eddington accreting 
stellar-mass black hole, however, the slop e is steeper than 
that predicted bv IPoutanen et all (|2007| ). 

The spectral fitting with a p-free model is consistent 
with the L oc T 4 relation. A recent Suzaku observation 
of NGC 1313 X-2 showing a p-free disk temperature of 
1.86 keV with a lumi nosity of 10 40 ergs s" 1 at a given 
distance of 4.13 Mpc (Miz uno et al.|[2007h . however, de- 
viates from the best-fit L oc T 4 line from XMM-Newton 
observations, which could be caused by the considerable 
error on the temperature or the difference between in- 
strumental calibrations. In the L — kT diagram (Fig. [4] 
bottom), NGC 1313 X-2 appears on the L oc T 4 exten- 
sion of the L—kT patterns of Galactic black hole binaries. 
Again, this favors an interpretation in which the source 
is a stellar-mass black hole undergoing super-Eddington 
accretion. Since the normalization of the L oc T 4 line is 
dependent on the black hole spin, color correction fac- 
tor, and temperature profile correction factor, the mass 
of the compact object cannot be precisely determined 
from Figure |U We estimate the compact object mass of 
NGC 1313 X-2 to be around 10-20 M Q as compared to 
Galactic black holes. 

We note that these results on NGC 1313 X-2 do not 
exclude the possibility that the soft excesses in other 
ULXs may be due to accretion disks. As described in 
iFeng &; Kaaretl (|2006f ) , the spectral evolution — even con- 
sidering only the relation of the photon index versus lu- 
minosity of NGC 1313 X-2 — is unlike that seen from stel- 
lar mass black holes. Repeated observations of several 
ULXs will be required to measure the spectral variabil- 
ity patterns of a sample of ULXs to determine if the 
behavior of the soft excess is, in general, consistent with 
that expected for accretion disk emission. 



We thank Marek Gierlinski for providing the data 
of Galactic black hole binaries and the anonymous ref- 
eree for the insightful comments, particularly regarding 
the uncertainty in the MCD temperature determination, 
which improved the paper significantly. 
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